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We carry out a self consistent calculation of the structure 
functions in the dissipation range using Navier Stokes equa- 
tion. Combining these results with the known structures in 
the inertial range, we actually propose crossover functions for 
the structure functions that takes one smoothly from the in- 
ertial to the dissipation regime. In the process the success of 
the extended self similarity is explicitly demonstrated. 

PACS : 47.27.Gs, 47.27.Eq, 05.45.+b, 05.70.Jk 

The inertial range of fully developed homogeneous, 
isotropic turbulece has been investigated extensively [1-8] 
in the past decade. In comparison far dissiaption range is 
been less well studied and as far as we know, a systematic 
study of the structure functions, based on Navier Stokes 
equation (NS), has not been carried out. In this work, we 
report a self consistent calculation of the structure func- 
tions in the dissipation range. Using accepted results in 
the inertial range, we propose crossover functions for the 
structure functions and thus demonstrate how extended 
self similarity can be understood. 

We work with forced three dimensional NS equation 
for incompressible flows, written in the momentum space 
as, 

^(k) + ^fc 2 w 4 (k) = yAf i ji(k)53*i J (p)»i(k-p) + / i (k,t) 

p 

(1) 

Where My/ = kjPa(\s.) + fc;Py(k) and the transverse 
projector, Pij = S%j — kikj/k 2 , where the external noise 
/, is S correlated and is necessary to maintain the en- 
ergy balance in the inertial range. The energy input per 
unit time (e) at the long wave lengths cascades through 
different lenght scales due to the nonlinear term and for 
k > kjj, is dissipated by molecular viscosity (u), here 
k D = (e/i/ 3 ) 1 / 4 . For (k « k D ), we have the so called 
inertial range, where one expects, 

S 2n ~ fc-^"+ 3 ") (2) 

with structure function defined as in Dhar et al jj| 
as, S n = (\v(k)\ n ). The exponent ( n is n/3 in the Kol- 
mogorov limit. In general it differs from n/3 and one 
of the best estimates of the deviation is due to She and 
Lcvcquc |l(J which gives, 



£C„ = Cn-rV3 = -^ + 2[l-(^)"/ 3 ] (3) 

In this work we investigated the dissipation range and 
our principal results are 

S 2n (k) ~ k n ^e- nk l K , (k»K) (4) 

where 5 2 = 2 — D (D being the dimensionality of space) 
and K = 0(/cd). By studying the correction to the above 
result as powers of K/k, we propose (in D = 3) the 
crossover function (crossover from dissipation to inertial 
range) 

S 2 (fc)~i(l + aq|) 2 +< 2 e- fc /* (5) 

where oq is number of 6(1), while for higher order 
structure function, 

S 2n (k) ~ k nS2 (l + a n )"(3+*0+&» e -n*/JC ( 6 ) 
k 

The constants a n are non-universal but will be shown 
to be almost independent of n. The explicit crossover 
forms that we have written down helps us understand 
the idea of extended self similarity (ESS) introduced by 
Benzi et al pj. Our approach is alternative to that of 
Segel et al jF3j. Writing Eq.6 by expanding about the 
inertial range form, we note that, 

S 2n (k) ~ fc-(C-+3n) (1 _ [n "O + ^+Can ^ (?) 

a n K 
the simple power law will break down when, 



n(a„ —3 — S2) — 

From phenomenology, it is known that S n falls off from 
the fc _ (C2n+3n) i me m dissipation range. This con- 
strains a n (from Eq 7 and using 5 2 — 2 — D), 

n(a n - 2) - C 2 » Q 

Now if we assume a n > and use the fact that a n is 
almost independent of n (shown later), we get, 

ai > 2 + C2 

As a n — 2 > C2(= 2/3), the difference n(a n — 2) — <^ 2n 
will grow with n (since we know that C, 2n deviates more 
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from linearity for higher moments). Hence for higher n 
the S n curves will fall off from the scaling regime at even 
lower k values. This is completely consistent with the 
standard phenomenology filjfl . 

We now turn to ESS. From Eq.6 it is clear that, 

S 2n (k) ~ [S 2m m [ t^ ] e- nk,K e m{ ^k )k/K 

[l + q-l*]"(3+fr)+C,„ 

X [l + a™ 1 #] (m(3+52)+C2m)( ^^ ) 

The explicit k dependant terms on the r.h.s. of the 
above expression will cause deviation from scaling. But 
it is apparent that the exponential factor is much more 
weakly decaying than e ~ nk / K (i n fact it is constant for 
the Kolmogorov situation of ( 2n = 2n/3) and also the 
variation of [1 + a" 1 JL] n ( 3 + s *)+b» j s mu ted by the de- 
nominator (as a„s' have been assumed to have the same 
sign and shown to be almost independent later). Con- 
sequently a plot of log^n vs log will remain 
a straight line over a far longer range than S 2n vs 
fc -(C2„+3n)_ Tnis is tne con tent of ESS. Few other phe- 
nomenological consequences are also manifest. For ex- 
ample, as (n — to) grows the scaling regime will become 
gradually shorter. In fact with a n independent of n to 
a first approximation and £2™ almost proportional to n, 
the scaling of log S 2n vs log S 2m is virtually exact. 

We first note that correlation functions in the dissipa- 
tion range falls off extremely fast jl^] with the charac- 
teristic scale ko and because of the existance of the scale 
there is no divergence in the self energies and correlation 
functions. Absence of divergence in the self energy im- 
plies that vicosity coefficient v is not renormalised. The 
correlation function is given at the self consistent single 
loop level by, 

S 2 (k, w) = \G\ 2 k 2 J ^^ a (k, p, k p) 

xS 2 (\k-p\,w-w')S 2 {p,w') (9) 
where the angular factor, 

a(k, p, k - p) = -(1 - xyz - 2y 2 z 2 ) 

The trio (k, p, k — p) forms a triangle and x, y, z 
are the direction cosines of the angles opposite to k, p 
and k — p respectively. The response function G _1 = 
—iw + vk 2 and the correlation function S 2 {k, w) — 
k 52 f(k/k D )k 2 /(w 2 + v 2 k A ) 7 such that 

/T 
^S 2 (k,w) = S 2 (k,t = 0) = k 5 *f(k/k D ) (10) 



Comparing the two sides of Eq.9, the function f(k/kn) 
has the structure f(k/ko) = e~^ k l kD , since on the right 
hand side of Eq.9, 

e -f3p/k De ~t3\k-p\/k D = e -/3(|!-q| + |#+q|)/fc D 

= e -0k/k D r e -e{q 2 /k 3 ) + N 

~ e - pk/kD (11) 

thus reproducing the exponential factor of the left hand 
side. Power counting of the momentum in Eq.9 now leads 
to, 

S 2 = -(£>- 2) (12) 

To check the correctness of this self consistent solu- 
tion we avaluated the equal time limit of the integral 
on the r.h.s. of eqn.9 numerically (with lower cutoff kp 
and D = 3). In fig. 1(a) we plot this integral Id(k) and 
compare it with the function -j-e~ fe /' CD . The agreement is 
good for k/kn > 20 (ie, k » ko)- 

If the above formalism has to approach the crossover 
behaviour, then we need to include the first correction 
to the large k behaviour. We do this by saying that the 
correction is in powers of (K/k) m and thus in D = 3 

S 2 = h -^[l + b 1 {K/k) m ]e- k l K (13) 

where K — fco//3. The right hand side of Eq.9, lin- 
earised in bi and considered at zero frequency can be 
written as being proportional to: 

d 3 p ^ e -( P +|k-p|)/jf 

J(^)3«(k, P ,k P) (p2 + | k _ pmk _ p | 

x[l + b 1 (K/p) m + h(K/\k - pi)" 1 ] (14) 

The requirement that the integral involving b\ is finite, 
leads to m < 2. If we write m = 2 — e, we can evaluate 
the integrals to the leading pole |Q in e. The integral 
not involving b\ can be evaluated using saddle point tech- 
nique, with dominant contribution coming from p ~ k. 
The result of the above manipulation must be of the form 
[1 + bi(K/p) m ] for self consistency and at the level of 
approximation just described, we find m = 1. Thus, we 
have the result that for k » k D , S 2 oc ^(l+bi^)e~ k ^ K 
and for k « kjj (inertial range), S2 oc fc - ( 3+ ^ 2 ). The 
simplest interpolation is Eq.5. 

Using the above analysis as a guide towards determin- 
ing m, we numerically evaluated the correction integral 
I c (k) of 9 (61) in eqn.14. Using the same values of bo 
and the lower cutoff kjj which we had used for fitting the 
dominant term, we find self consistency of the correction 
integral can be achieved for to ~ 1/4. In fig. 1(b) we plot 
this integral I c (k) as a function of k/ko and compare it 
with the 0(6i) term in eqn.13. The agreement is good 
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for k/ko > 7. We have chosen K — kp for our numerics. 
All our arguments demonstrating the ESS properties of 
S„ hold good as long as m > 0. 

It should be noted that in this far dissipation range 
that we are considering here, the single loop self consis- 
tency is sufficient. We have checked that the contribu- 
tions from higher (> 2) loop diagrams are at most of the 
same order as the single loop diagram. So their inclusion 
just changes the amplitude of 82(h). This statement is 
true for the evaluation of S^n ( n > 1) also, which we do 
now. Out of the various possible arrangements of the k 
and — k external legs on an one loop diagram, we evalu- 
ate the most relevant one (shown in Fig.l). Contribution 
from other possible one loop diagrams are exponentially 
smaller and hence their contributions are negligible. The 
contribution from Fig.l is, 

■MM) = <KMM-M)]"} 







k 2n 1 dh.. 


/ dt 2n \ 


J — OO 


J — OO J 



xG(k;t,ti)G(-k;t,t 2 ) ■ ■ ■ 

xS 2 ( P , |*i - t 2 |)5 2 (|k - p|, \t 2 - t 3 \) . . . (15) 

As S 2 (M) ~ fe-^-^e-^e-l*!^ 2 , we note that the 
integral of Eq.15 will be dominated by the low momen- 
tum pole at p — k. Using a pole approximation for eval- 
uating the integral, a momentum count produces the re- 
sult that S 2n (k,t) oc k n&2 e~ nk / K . This establishes Eq 
4. However, within this formalism, though we cannot 
rigorously show that Eq.4 holds for odd moments also, 
for monotonicity sake we assume this to be true. Now 
we note that Eq 4 implies S n ~ ie, simple scaling 
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FIG. 2. Diagram contributing to S2n(k) in the dissipation 
regime 

behaviour results in the far dissipation range. This is in 
mild contrast to the simulation results ||, where very 
weak multiscaling (ie, very small deviation from n/3) 
has been reported. But given that this deviation is very 
small, e.g. for n = 7 the numerical exponent is 2.24 in- 
stead of 7/3, our estimate for this far dissipation range 
is a very close one. 

We now study the first deviation of 5*2n from its form 
in Eq 4. To do so, we introduce the first deviation of 
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S2n{k,t) in Eq 15. The integral in in Eq.15 is already 
pole dominated and hence the additional part is pole 
dominated as well. There are n contributions of equal 
strength from each of the Ssip) and S^Gk — p|) and con- 
sequently for k >> K 



S 2n oc k- nl ^ D - 2 \\ + b n (K/k)]e- nk / K 



(16) 



where b n oc nb\. With the quantity n(3 + 62) + C2n 
in Eq.6 roughly proprtional to n, we consequently infer 
that in the interpolation formula of Eq.6, the constant 
a n is to a good approximation independent of n. Thus 
the main results Eq.4 - Eq6, are obtained. 

Now we look at the real space structure function 
S^M = ([ v ( x + r ) — v ( x )] 2 ) which is the inverse Fourier 
transform of 2[uQ 2 5(k) — £2 (A;)] (where uo 2 /2 is the mean 
energy). For r in the far dissipation range S^r) will be 
determined by our k >> form of S^fc) (ie,~ ^e~ k l K ). 
This yields S2 (r) = c\r 2 + 0(r 4 ). Here c\ is a function 
of v, e. This form of S 2 (r) is consistent with the result 
of Sirovich et.al. |16|. The added advantage of our k— 
space calculation is the ability to predict the higher order 
structure functions (S l 2n(fc)) also. 
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In summary we have shown that by considering Navier 
Stokes equation and doing a self consistent treatment 
of the dissipation range (characterised by the existence 
of a scale krj ), we can establish forms for the vari- 
ous order structure functions. By the first correction to 
the asymtotic situation and using the known results in 
the inertial range (k << kjj), we can construct explicit 
crossover functions for the structure functions (crossover 
from k » kr> to k « kr> ) . The validity of ESS is easy 
to see. 
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